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Abstract 

A three-dimensional stress-concentration analysis 
was conducted on straight-shank and countersunk 
(rivet) holes in a large plate subjected to various 
loading conditions. Three-dimensional finite-element 
analyses were performed with 20-node isoparametric 
elements. The plate material was assumed to be 
linear elastic and isotropic, with a Poisson’s ratio 
of 0.3. Stress concentrations along the bore of the 
hole were computed for several ratios of hole radius to 
plate thickness (0.1 to 2.5) and ratios of countersink 
depth to plate thickness (0.25 to 1). The countersink 
angle was varied from 80° to 100° in some typical 
cases, but the angle was held constant at 100° for 
most cases. For straight-shank holes, three types 
of loading were considered: remote tension, remote 
bending, and wedge loading in the hole. Results 
for remote tension and wedge loading were used to 
estimate stress concentrations for simulated rivet or 
pin loading. For countersunk holes, only remote 
tension and bending were considered. Based on the 
finite-element results, stress-concentration equations 
were developed. Whenever possible, the present 
results were compared with other numerical solutions 
and experimental results from the literature. 

Introduction 

A riveted joint is a commonly used method of 
joining structural components. Joining introduces 
discontinuities (stress risers) in the form of holes, 
changes in the load path due to lapping, and ad- 
ditional loads such as rivet bearing and bending 
moments. Because of these changes at the joint, 
local stresses are elevated in the structural compo- 
nent. Accurate estimations of these local stresses are 
needed to predict joint strength and fatigue life. 

Exhaustive studies on stress-concentration factors 
(SCF’s) for holes and notches in two-dimensional 
(2-D) bodies subjected to a wide variety of loadings 
have been reported in the literature (refs. 1 and 2). 
Studies have also been made on three-dimensional 
(3-D) stress concentrations at circular holes in plates 
subjected to remote tension loads (refs. 3 to 6). A 
recent paper by Folias and Wang (ref. 6) provides a 
review of these previous solutions and presents a new 
series solution. The Folias and Wang solution covers 
a wide range of ratios of hole radius to plate thick- 
ness. The stress concentration at a hole in a plate 
subjected to bending was first presented by Neuber 
(ref. 4) using the Love- Kir chhoff thin-plate theory 
(ref. 7). Reissner (ref. 8) rederived the plate solution 
including the effect of shear deformation and showed 
that Neuber’s solution was unconservative. Reiss- 
ner’s SCF solution for bending loads is presented in 


terms of Bessel functions. Naghdi (ref. 9) extended 
Reissner 7 s analysis to elliptical holes using Math- 
ieu’s functions. Rubayi and Sosropartono (ref. 10) 
conducted 3-D photoelastic measurements to ver- 
ify Reissner’s circular hole and Naghdi’s elliptical 
hole solutions. Many other analytical (e.g., refs. 11 
and 12) and experimental (e.g., refs. 13 to 15) results 
are reported in the literature for remote loading, but 
none consider 3-D effects for rivet loading in the hole. 

Only two papers in the literature report results 
on stress concentration at countersunk holes (refs. 16 
and 17). Both papers use the 3-D photoelastic slice 
method to obtain stress concentrations for holes in 
thick plates. Cheng’s results (ref. 17) include stress- 
concentration factors for both tension and bending 
loads. 

Three-dimensional stress concentrations at riv- 
eted joints are not fully understood. Knowledge of 
3-D stress concentrations is needed to verify the ad- 
equacy of 2-D solutions in the handbooks. Further- 
more, 3-D stress concentrations are needed to pre- 
dict the strength and life of joints. Countersunk-rivet, 
construction is commonly used in aircraft industries 
to achieve acrodynamically smooth surfaces. These 
joints are not amenable to 2-D approximation, and 
stress concentrations for countersunk holes have not 
been reported in handbooks. 

The objective of the present study is to conduct 
a comprehensive analysis of three-dimensional stress 
concentrations for circular straight-shank and coun- 
tersunk (rivet) holes in a large plate subjected to 
various loads encountered in structural joints. Three 
types of loading, remote tension, remote bending, 
and wedge loading in the hole (simulated pin load- 
ing), are considered for the straight-shank hole. Two 
types of loading, remote tension and remote bending, 
are considered for countersunk holes. 

Three-dimensional finite-element (F-E) stress 
analyses of large plates with straight-shank and 
countersunk circular holes were conducted with the 
FRAC3D F-E code. The FRAC3D code is an elas- 
tic and elastic-plastic fracture mechanics code de- 
veloped at NASA Langley Research Center for the 
analysis of cracked isotropic or anisotropic solids, 
based on the 20-node isoparametric element. A wide 
range of hole sizes (ratio of hole radius to plate 
thickness) and countersink depths (ratio of coun- 
tersink depth to plate thickness) are considered in 
the analysis. The plate material is assumed to be 
linear elastic and isotropic, with a Poisson’s ratio 
of 0.3. The influence of the countersink angle on 
stress concentrations is also examined. With the 
F-E results, simple, series-type stress concentration 



equations are developed for a wide range of hole- 
radius-to-plate- thickness ratios and for any location 
along the bore of the hole. Whenever possible, the 
present results are compared with results from the 
literature. 

Nomenclature 


b 

depth of straight-shank portion of hole 

FEM 

finite-element method 

h 

one-half height of plate 

Kb 

stress- concentration factor along bore 
of hole under bending 

Kb 

(Amix 

maximum stress-concentration factor 
along bore of hole under bending 

K v 

stress-concentration factor along bore 
of hole under pin loading 

K, 

stress-concentration factor along bore 
of hole under tension 

^niax 

maximum stress-concentration factor 
along bore of hole under tension 

K w 

stress-concentration factor along bore 
of hole under wedge loading 

M 

applied remote bending moment 

P 

applied pin or wedge loading 

r 

radius of straight,-, shank portion of 
hole 

S 

applied remote tension stress 

t 

plate thickness 

w 

one-half width of plate 

x,y, z 

Cartesian coordinate system 

a ij 

coefficients in stress-concentration 
equations 

ft; 

coefficients in stress-concentration 
equations 

9c 

countersink angle 

a yy 

hoop stress at </> — 90° 

<P 

angle defining applied stress distribu- 
tion in hole 


Rivet Hole Configurations 

Two types of rivet hole configurations, straight- 
shank and countersunk holes in a large plate, were 
considered. The two types of holes and nomenclature 
used are shown in figure T In the plate with a 

2 


countersunk hole, the thickness was divided into 
two sections: the cylindrical section, referred to as 
the straight-shank depth 6, and the conical section, 
referred to as the countersink depth t-b. The two 
sections meet to form an edge referred to as the 
countersink edge. The stress concentrations depend 
on the length of the straight-shank or countersink 
depth. The two extreme cases of countersunk holes 
are when b - 0 (knife edge) and b = t (straight-shank 
hole). The 1 countersink angle was 9 C . (See fig. 1(b).) 

In the present study, the plate width and height 
were selected large enough so that the stress- 
concentration solutions were not greatly affected by 
the remote boundaries. Stress concentrations for 
finite-size plates have to be generated with the use of 
cither analytical or numerical methods. A wide range 
of values for h ole- radi us- to-plate- thickness ratio r/t 
and straight.-shank-depth-to-plate thickness ratio b/t 
were considered in generating the data base on stress- 
concent ration solutions. For the straight-shank hole, 
six values of r/t (encompassing the range of struc- 
tural configurations used in industry) were selected: 
2.5, 1.5. 1.0, 0.5, 0.25, and 0.1, with w/r — 5 and 
h/r = 5. For the countersunk-hole configuration, r/t 
values selected were 2.0, 1.0, 0.5, and 0.25 and b/t 
values were 0. 0.25, 0.50, and 0.75 with w/r = 7.5 
and h/r -- 7.5. Although the straight-shank hole 
configuration is a special case of the countersunk hole* 
(b/t — 1), for convenience the two configurations are 
considered l separately. 

Loading Conditions 

Figure 2(a) shows the three types of loadings 
that were applied to the plate with the straight- 
shank hole?: remote tension stress S , remote bend- 
ing moment, per unit width A/, and wedge load- 
ing P. Appendix A explains how the stress 
concentrations for a pin-loaded hole were approxi- 
mated from the remote tension and wedge loading 
solutions. The wedge loading was imposed on the 
hole boundary as a normal pressure loading that has 
a cosine distribution and is assumed to be constant 
through the plate thickness. The surface pressure is 
defined as (2P/wrt)cos 4> (refs. 18 and 19) and was 
applied over the angle <p — ±90°. The angle & is 
measured from the y-axis. (See fig. 2(a).) The bend- 
ing moment. M was applied as an equivalent remote 
stress that, varies linearly through the plate thick- 
ness. For countersunk holes, two loading types, re- 
mote tension and remote bending, were considered. 
(See fig. 2(b).) Because of the lack of understanding 
of 3-D load transfer between the rivet and the coun- 
tersunk hole, rivet (pin) loading was not considered 
in the current study. 
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Definition of Stress-Concentration Factor Finite-Element Modeling 


Although the definition of the stress- concent rati on 
factor is given in many classical books on theory 
of elasticity and in stress-concentration handbooks, 
many of these solutions are associated with 2-D con- 
figurations. For 3-D configurations, however, the 
stress concentration varies along the structural dis- 
continuity, such as along the bore of the hole. Herein, 
the stress-concentration factor is defined as the stress 
at any location along the bore of the hole normalized 
by a characteristic stress (related to applied loading). 
For configurations and loading conditions considered 
in this study, the highest stresses occurred along the 
bore of the hole at the intersection of the hole surface 
and the y — 0 plane. Even for the case of pin loading, 
the peak stresses occurred at (f> — 90° because the pin 
contact angle was assumed to be 90°. (See appen- 
dix A for details.) The stress-concentration factors 
for the three loading conditions arc defined as follows. 

Remote tension . The stress-concentration fac- 
tor for tension Kf is the hoop stress a yy at (j> = 90° 
along the bore of the hole normalized by the applied 
remote tension stress 5 and is given by 

KM = "Jp ( 1 ) 


Remote bending . The stress-concentration fac- 
tor for bending K ^ is the hoop stress (Tyy at <fi = 90° 
along the bore of the hole normalized by the remote 
outer-fiber bending stress 6A//£ 2 and is given by 


Kb(*) 


°yy(z) 

6 M/t 2 


(2) 


Wedge loading . The stress-concentration factor 
for wedge loading K w is the hoop stress a yy at 
cj) — 90° along the bore of the hole normalized by 
the average bearing stress P/2rt and is given by 


K w (z) 


^yyjz) 

P/2rt 


(3) 


Pin loading . The stress-concentration factor for 
pin loading K p is obtained from a superposition of 
remote tension and wedge loading. (See appendix A.) 
The factor K p is defined as the hoop stress <r yy at 
<f> = 90° along the bore of the hole normalized by the 
average bearing stress P/2rt and is given by 


K P (z) 


P/2rt 


(4) 


A three-dimensional finite-element, code FRAC3D 
developed at NASA Langley Research Center for an- 
alyzing cracked isotropic and anisotropic solids was 
used in this study. The code is based on the 20-node 
isoparametric element formulation. The stiffness ma- 
trix and the consistent load vectors were generated 
with the 2 by 2 by 2 Gaussian quadrature formula. 
The program uses a vector skyline Choleski decom- 
position algorithm (ref. 20) for solving matrix equa- 
tions of equilibrium. The plates with the straight- 
shank hole and remote tension and wedge loading 
were symmetric about the x — 0, y = 0, and z = 0 
planes. The remote bending was symmetric about 
the x — 0 and y — 0 planes and antisymmetric about 
the 2 — 0 plane. Because of these conditions, only 
one-eighth of the straight-shank hole plate was mod- 
eled. The FRAC3D code has an option to impose 
symmetry and antisymmetry boundary conditions. 
The plate with the countersunk hole was symmetric 
about the x = 0 and y = 0 planes; hence, one-fourth 
of the plate was modeled. The F-E model includes 
the full thickness of the plate. 

Because many configurations were to be analyzed, 
a simple 3-D modeling procedure was developed to 
generate the finite-element, meshes. In this proce- 
dure, a 2-D F-E mesh in the x-y plane was gener- 
ated with refined elements near the hole boundary. 
Then the 2-D mesh was translated in the z-direction 
(with appropriate x-y transformation to account for 
the countersunk hole). Typical 3-D F-E meshes for 
one-eighth of a straight-shank hole in a plate and 
for one-quarter of a countersunk hole in a plate are 
shown in figure 3. 

For all straight-shank hole models, the half- 
thickness of the plate was divided into six layers of 
unequal thickness. The layer thicknesses (starting 
from the z — 0 midplane) were 15, 13, 10, 6, 4, 
and 2 percent of the total plate thickness. The small 
thickness layers were used in the high-stress-gradient 
regions (near the free surface). The F-E model had 
936 elements and 4725 nodes (14 175 degrees of free- 
dom). For different values of r/i, the hole radius 
was kept constant and the plate thickness was scaled 
by t/r. The F-E mesh for rjt — 1.0 is shown in 
figure 3(a). 

In the countersunk hole, there are three regions 
where the stress gradient is high: near the two free 
surfaces of the plate and at the countersink edge. 
Therefore, different through-the-thickness idealiza- 
tions were used for different countersink edge loca- 
tions b/t. Table 1 gives the details of the F-E ide- 
alizations used for countersunk holes with b/t — 0, 
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0.25, 0.50, and 0.75. Figure 3(b) shows a typical F-E 
model (r/t = 0.25 and b/t = 0.50) for one-quarter of 
a plate with a countersunk hole. 

Comparison With Other Solutions 

The present 3-D stress-concentration factors 
(SCF’s) for the straight-shank hole arc compared 
with Folias and Wang’s solution (ref. 6) for re- 
mote tension and with Reissner’s solution (ref. 8) 
for remote bending. Three-dimensional stress- 
concentration solutions for wedge loading or simu- 
lated pin loading have not been reported in the litera- 
ture. For countersunk holes, the present solutions are 
compared with Cheng's photoelastic measurements 
(ref. 17) for thick plates subjected to tension and 
bending. 

Straight-Shank Hole 

Remote tension . The distribution of the stress- 
concentration factor Kf along the bore of the hole for 
remote tension is shown in figure 4 for various values 
of r/t. The stress concentrations are symmetric 
about the inidplane (z/t = 0). (Note the expanded 
scale on the ordinate axis.) In all these cases, the 
plate width and height were selected large enough 
( w/r = hjr = 5) so that k't values are not greatly 
affected by the finite plate. The SCF’s for r/t < 0.5 
are about 2 percent larger than Folias and Wang’s 
(ref. 6) infinite-plate solutions for all 2 -values. (For 
clarity, results from ref. 6 are not shown in fig. 4.) 
Part of this difference may have’ been caused by 
the finite-size plate used in the present study. For 
r/t = 0.1, the present results show the same trend 
as that of the Folias and Wang solution (maximum 
SCF near the free surface), but the magnitude of the 
present results is 3 percent lower than the magnitude 
of their solution. It is expected that the classical 
value of Kf = 3 would be obtained for much thicker 
or thinner plates. From the plots of Kt for various 
values of r/t, two observations are made: 

1. For r/t > 0.5, the maximum K% occurs at 
z/t, — 0 (midplane). For thicker plates ( r/t < 0.25), 
the maximum Kt location shifts toward the free 
surface (z/t = ±0.5). 

2. At z/t = 0, the Kt value appears to peak at a 
value of 3.22 for r/t = 0.5. 

These two trends are consistent with crack-tip stress- 
intensity factors for cracks in thick plates, as ob- 
served by many investigators. (See, for example, 
ref. 21.) However, such a drop in Kt near the 
midplane (z/t = 0) for very thick plates was not 


reported in reference 6, even for r/t as small as 0.02. 
The reason i'or not capturing this expected trend in 
reference G is unknown. 

Remote bending . Figure 5 shows a compar- 
ison of maximum bending stress-concentration fac- 
tor Kh calculated from the present analysis, from 
Reissner’s shear deformation plate theory (ref. 8), 
and from Neuber’s thin- plate theory (ref. 4). Because 
both Reissner and Neuber assumed that the stress 
distribution was linear through the thickness, K\ } 
always occurs at z/t = ±0.5. However, the present 
F-E solutions show that the location of Kh is at 

fJ rnax 

z/t — ±0.5 (free surface) for thin plates (r/t > 0.5), 
but the maximum SCF is slightly interior from the 
free surface’ (\z/t\ < 0.5) for thick plates (r/t < 0.25). 
(These results are shown subsequently.) In figure 5. 
the maximum SCF values from the F-E analysis 
arc plotted for various values of r/t extrapolated 
to r/t — 0. (See the dashed curve.) Results for 
Reissner’s shear deformation theory and the present 
results agree well with each other for r/t > 1.5. 
The difference between results for Reissner’s solu- 
tion and the present results for r/t less than unity 
is about 4 to 8 percent. Neuber’s thin- plate the- 
ory. Kf, — (5 ± n)/(3 ± j /), is inadequate even for 
r/t — 2.5 and produces values about 6 percent lower 
than those for Reissner’s solution and the present 
results. 

Countersunk Hole 

Chong (ref. 17) measured 3-D stress-concentration 
factors for countersunk holes in thick plates using 
a photoelastic slice technique for both tension and 
bending. Cheng’s photoelastic models for tension 
(model 7) and for bending (model 8) were analyzed 
through the generation of separate F-E meshes. The 
geometric parameters of models 7 and 8 are given in 
table 2. For both models, b/t = 0.6 and 0 C — 90°. 
The SCF’s for the two configurations at the critical 
locations are presented in table 2. The F-E results 
show that the maximum SCF for remote tension oc- 
curs slightly away from the countersink edge and in 
the straight-shank portion of the hole (at z/t = 0.08, 
whereas the countersink edge is at z/t — 0.1). The 
maximum SCF calculated from the F-E analysis is 
within 3 percent of Cheng’s measured value. (See 
table 2: note that percent error is defined as the 
difference between solutions divided by the largest 
stress-concentration value.) For bending, three lo- 
cations on the hole (z/t — —0.5, 0.1, and 0.5) were 
considered for comparison. The difference between 
Cheng’s measurements and the present solution is 
about 2.5 percent at z/t = —0.5, but the difference 



is about 8 percent at the countersink edge ( z/t = 0.1) 
and along the countersink flank (z/t = 0.5). As 
previously observed for straight-shank holes in thick 
plates, the maximum SCF is not at the free surface 
(z/t = —0.5) but is slightly interior to the free surface 
(z/t = —0.48). The drop in SCF at the free surface 
is attributed to the well-known free- bound ary- layer 
effect (refs. 21 and 22). 

Effect of Countersink Parameters 
on SCF 

The two parameters that can influence the SCF 
for countersunk rivet holes are the countersink angle 
0 r and the countersink depth t — b. (See fig. 1(b).) 
The effects of these two parameters on SCF at coun- 
tersunk holes in plates subjected to tension and bend- 
ing were analyzed. 

Countersink Angle 

The effect of small variations in the countersink 
angle 0 r on tension and bending SCF was analyzed 
with Cheng’s model 7 configuration (ref. 17). Fig- 
ure 6 shows the distribution of Kt and K ^ along the 
bore of the hole for 0 C = 80°, 90°, and 100°. A change 
in 0 ( of ±10° from the reference angle of 90° changes 
maximum K\ by about 3.5 percent at the countersink 
edge. However, the variation in Kt is much smaller 
at all other locations on the hole boundary. For the 
±10° variation in 0 C , K\ } varies less than 1 percent. 
(See fig. 6(b).) These results are for a thick plate, 
where r/t — 0.24. For thin plates, used in aircraft 
applications (r/t of about 2), the effect of 0 C varia- 
tion on SCF is of the same order as that shown for 
the thick plates. 

Countersink Depth 

As previously mentioned, the countersink-depth- 
to-plate-thickness ratio is defined as 1 — (b/t), where 
b/t represents the ratio of the straight-shank depth 
to the plate thickness. For convenience, b/t is used as 
a depth parameter. Figure 7 shows the distribution 
of Kt and K\> along the bore of the hole (—0.5 < 
z/t < 0.5) for b/t = 0, 0.25, 0.50, and 0.75. The 
plate w/r = h/r = 7.5, r/t = 2.0, and d c — 100° (a 
typical value for the aircraft industry). In figure 7(a), 
the maximum tension SCF’s occur at the countersink 
edge for all values of b/t. However, in thick plates, 
the maximum SCF’s occur slightly away from the 
countersink edge, on the straight-shank portion of 
the hole. The maximum tension SCF’s are 4.06, 
4.10, 3.82, and 3.39 for b/t = 0, 0.25, 0.50, and 0.75, 
respectively. The highest SCF for the countersunk 
hole under remote tension is about 37 percent higher 
than the classical 2-D value (Kt = 3). In contrast, for 


finite-thickness plates, SCF’s for countersunk holes 
are only about 30 percent higher because the finite 
thickness elevates the SCF, as shown in figure 4. 
Figure 7(b) shows the bending SCF distribution for 
various values of b/t. In contrast to tension loading, 
the maximum bending SCF is almost unaffected by 
b/t , except for b/t — 0 (near the knife-edge location). 
The variation of maximum K{> (at z/t — ±0.5) is less 
than 1 percent for 0.25 < b/t < 0.75. The maximum 
bending SCF at z/t. — —0.5 is —2.41 for b/t = 0; this 
Kfr is about 24 percent lower than that for b/t = 0.50. 
(Note that the SCF will be positive at z/t = —0.5 if 
the moment is reversed.) 

Stress-Concentration Factor Equations 

In this study, 3-D stress-concentration factors 
for a wide range of hole configurations and load- 
ings were generated with the finite-element method. 
These solutions may be used in structural design 
as they are or they may be interpolated to cal- 
culate stress concentration at any other location 
along the hole boundary or for other hole (r/t and 
b/t) configurations. The F-E solutions will be eas- 
ier to use if equations are developed. In this sec- 
tion, SCF equations are developed by fitting the 
F-E results to double-series polynomial equations. 
Separate equations are developed for straight-shank 
and countersunk holes subjected to different loading 
conditions. 

Multiparameter least-squares equation fits were 
performed with the International Mathematical and 
Statistical Library routine QRASOS, which uses the 
Householder transformation for matrix factorization 
(ref. 23). The weight factor for each of the SCF 
values along the bore of the hole is selected such that 
the weight is proportional to the length between the 
two neighboring points on either side of the point 
under consideration. For example, a weight factor 
for the ith point is (z r +\ — z/) /2 1. This procedure of 
selecting the weight factor minimizes the area under 
the SCF curve on the z- axis. Also, this procedure 
gives good fits to even unequally spaced data points 
without higher order oscillations, which are generally 
present in a high-order polynomial fit. 

Straight-Shank Hole 

The configuration for a plate with a straight- 
shank hole is symmetric about the 2 = 0 plane. 
(See fig. 2(a).) The tension and the wedge loading 
are symmetric about the 2 = 0 plane, whereas the 
bending is antisymmetric about the 2 = 0 plane. 
Therefore, an even-power polynomial in 2 and a 
general polynomial in r/t were used to fit SCF results 
for tension (Kt) and for wedge loading (K w ). An 
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odd-power polynomial in z and a general polynomial 
in rjt were used to fit the SCF results for bending 
(Kb). The forms of the SCF equations are 

4 4 

K <n = LE«ij( r /0 i W0 2j (5) 

t=0j=0 

where rn = t for remote tension and rn — w for wedge 
loading and 

4 4 

^ = ( 6 ) 

for remote bending. Equations (5) and (6) apply 
over the range 0.1 < rjt < 2.5. The orders of the 
polynomials for r/t and z/t are selected by trial and 
error such that the maximum difference between the 
F-E results and the equation results is minimal and 
the sum of residuals is a minimum. The coefficients 
o / j are given in tables 3, 4, and 5 for remote tension, 
wedge loading, and remote bending, respectively. 

Figures 8 to 10 show a comparison between the 
F-E results and the equation results for various val- 
ues of r/t for remote tension, wedge loading, and re- 
mote bending, respectively. (Note that an enlarged 
scale is used on the ordinate axis in figs. 8 to 10 
to magnify the difference between the equation and 
the F-E results.) The equation results agree well 
with the F-E results for all values of rjt and for 
all three loading conditions. The maximum differ- 
ence between the F-E solutions and the equation re- 
sults is about 1 percent. For both remote tension 
(fig. 8) and wedge loading (fig. 9), the SCF drops 
near the free surface. The drop is larger for thicker 
plates (smaller rjt). The bending SCF (fig. 9) is 
almost linear for rjt > 1.5 and becomes nonlinear 
for thick plates ( rjt < 1.0), particularly near the 
free surface. Thus, the assumption of linear stress 
(or strain) distribution through the thickness made 
in the Reissner (ref. 8) and Neuber (ref. 4) anal- 
yses is valid only for rjt > 1.0. For rjt. = 0.5 
and 0.25, the maximum SCF 1 is not at the free surface 
(z/t = ±0.5); it is located in the interior of the plate 
(I z/t\ < 0.5). 

Now that the SCF equations for remote tension 
and wedge loading have been established, the SCF 
equation for simulated pin loading is written as 

k „ = (7) 

Equation (7) is restricted to rjw = 0.2 because Kf 
and K w are generated for a plate with rjw — 0.2. The 


developmenl of equation (7) is given in appendix A. 
The results from equation (7) are shown in figure 11. 
Of course, these results show the same trends as 
those shown in figures 8 and 9 for tension and wedge 
1 o ad i ng . respect i vely . 

Countersunk Holes 

The configurations of the countersunk hole dic- 
tate that two separate SCF equations be fit: one 
equation for the straight-shank part (—0.5 < z/t < 

( b/t — 0.5)) and the other equation for the counter- 
sunk portion ((b/t — 0.5) < z/t < 0.5). Furthermore, 
separate^ equations were developed for each value of 
b/t. A general polynomial series equation in terms of 
rjt and z/t was fit to the F-E results with the 1 cast- 
squares procedure previously discussed. The SCF 
expiations are given by 

3 1 

Km = ± ± a ij( r / t Y (*/*■¥ ( 8 ) 

i= 0 j=0 

for —0.5 < z/t < (b/t — 0.5) and 
3 4 

K m - H W')' il z ~ b + (*/ 2 )] /(* - 6 )} J 

' 0 7-0 

(9) 

for (b/t) — 0.5 < z/t < 0.5. Equations (8) and (9) 
apply over 0.25 < r/t < 2.5. Coefficients a t] and (ijj 
for various values of b/t are given in tables 6 and 7 
for remote tension and remote bending, respectively. 
Figures 12 and 13 show comparisons between results 
for equations (8) and (9) with the F-E results for 
remote tension arid remote bending, respectively. 
The equation results and F-E results agree w r ell, 
except near the free surface for thick plates. Even for 
thick plates, the maximum SCF is wfithin 2 percent 
of the F-E results for all b/t values. Note that the 
bending SCF at z/t — 0.5 for the straight-shank hole 
(b/t -- 1.00. see fig. 10) is slightly less than that at 
z/t = —0.5 for the countersunk hole with b/t = 0.50 
(see fig. 13(c)). 

In appendix B, a FORTRAN program is given 
to evaluate the SCF’s for straight-shank and coun- 
tersunk holes subjected to remote tension, remote 
bending, pin loading, and wedge loading. This pro- 
gram is bast'd on equations (5) to (9) w r ith the coeffi- 
cients presented in tables 3 to 7. This program may 
be used to generate three-dimensional SCFs for any 
value of b/t and r/t and at any location along the 
bore of the hole. To generate the SCF's for values of 
b/t other than those used in this study, an interpola- 
tion scheme* between the available solutions has been 
implemented in the program. 



Concluding Remarks 

A comprehensive three-dimensional stress- 
concentration analysis of straight-shank and coun- 
tersunk (rivet) holes in a large plate subjected 
to various loading conditions encountered in ser- 
vice was conducted. The plate material was as- 
sumed to be isotropic, with a Poisson’s ratio of 0.3. 
Three-dimensional finite- element analyses were per- 
formed with 20-node isoparametric elements. Stress- 
concentration factors for wide ranges of hole-radius- 
to-plate thickness and countersink-depth-to-plate 
thickness ratio were generated. The countersink an- 
gle was varied from 80° to 100° in some typical cases, 
but tfie angle was held constant at 100° for most 
cases. For straight-shank holes, three types of load- 
ing, remote tension, remote bending, and wedge load- 
ing, were considered; for the countersunk hole only 
remote tension and remote bending were considered. 
Series- type equations were fit to the finite-element re- 


sults. These equations generally agreed within 1 per- 
cent of the finite-element results. 

Tension stress-concentration factor (SCF) for 
a countersunk hole was about 37 percent higher 
than the classical (2-D) solution for a circular hole 
(SCF = 3); the SCF was about 30 percent higher 
than the 3-D SCF for a straight-shank hole with the 
same hole-radius-to-platc-thickness ratio. However, 
the bending SCF was almost unaffected by counter- 
sinking the hole, except for the knife-edge case (no 
straight shank). Variation in the countersink angle 
(80° to 100°) had little effect on the peak SCF (a 
change of less than 3.5 percent) for both remote ten- 
sion and remote bending. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
April 15, 1992 
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Appendix A 

Computation of Stress-Concentration 
Factors for a Pin-Loaded Hole in a 
Large Plate 

This appendix shows how the SCF solutions for 
the wedge loading and for the remote uniform stresses 
can be used to predict the SCF for a pin-loaded 
(or rivet-loaded) plate. In this analysis, as already 
mentioned in the text, the pin and the plate surface 
are assumed to be smooth and the pin fits snugly into 
the hole (no clearance). Two-dimensional studies 
(refs. 24 to 27) have clearly demonstrated that the 
maximum tensile stress concentration due to pin 
loading occurs at the end of the contact between 
the pin and the hole boundary. For a snugly fit 
pin joint and modeling the pin, these references show 
that the contact angle is about 83°. However, Crews 
et al. (ref. 25) show that a contact angle of 83° or 
90° (between the pin and the hole) has very little 
effect on the maximum stress concentration (about 
a 3-percent difference). Because the present analysis 
assumes that the contact angle is 90°, the maximum 
stress concentration occurs also at 90°. Therefore, 
the superposition of the wedge loading and remote 
loading solutions (both give the maximum stresses 
at 90°) give the highest stress concentration. If 
the contact angle is assumed to be 83°, then the 
maximum stress at 90° is within 3 percent of the 
maximum value at 83°. 

In the present analysis, the pin-load reaction is 
approximated by a cosine distribution over the con- 
tact angle 0 of ±90° (refs. 18 and 19). (The angle is 
measured from the ?/-axis; see fig. 2(a).) Consider a 


pin-loaded hole as in the plate shown in figure 14(a), 
where the plate height h is large compared with the 
hole radius. This condition results in a uniform stress 
of P/2wt at y = —h. The SCF for this problem is 
measured by the (7 yy stress at y — 0 and x — r. Con- 
sider another problem, shown in figure 14(b), where 
the pin load acts on the lower half of the hole and the 
corresponding remote stress is P/2wt at y — h. The 
two problems in figures 14(a) and 14(b) are identi- 
cal except that the stress is 180° out of phase. On 
the x-axis, the a yy stress for these two problems is 
identical. The sum of these two loading conditions in 
figures 14(a) and 14(b) can be represented as wedge 
loading and remote tension, as shown in figures 14(c) 
and 14(d). Therefore, cr yy stress on the x-axis for the 
pin- load case is half the sum of the a yy stress due to 
wedge loading and remote tension. Alternatively, the 
stress concentration for the pin load K p is defined in 
terms of the stress-concentration factor for the wedge 
load K w and the remote tension load Kt as follows: 

Kr= K. + Ww) K L (A1) 

A plane strain analysis of r/w = 0.2 was performed 
with the 3-1) finite-element analysis described in 
the text. The computed pin-load SCF from equa- 
tion (Al) is K p — 0.994. This solution agrees rea- 
sonably well with the experimental results (0.985) 
reported by Chang et al. (ref. 24). Some differences 
are observed between the present results and those 
of other analyses, such as Crews et al. (0.87, ref. 25). 
De Jong (1.058, ref. 26), and Eshwar et al. (0.922, 
ref. 27). 
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Appendix B 

Computer Code Used To Calculate Stress-Concentration Factors 

This appendix presents a FORTRAN program and subroutine (SCF3D) used to calculate the stress- 
concentration factor at any location along a straight-shank or countersunk hole subjected to various loadings. 
This program was developed such that it may be readily incorporated into other stress-analysis or life-prediction 
codes. The program returns Kt, K^, K p , and K w for remote tension, remote bending, pin loading, and wedge 
loading. The simulated pin-load value K p was calculated from Kt , K w , and procedures that are presented in 
appendix A. 

PROGRAM MAIN 
C 

C COMPUTES THREE-DIMENSIONAL STRESS-CONCENTRATION FACTORS FOR 
C STRAIGHT-SHANK OR COUNTERSUNK HOLES SUBJECTED TO REMOTE 

C TENSION, REMOTE BENDING, PIN LOADING AND WEDGE LOADING 
C 

CHARACTER* 1 LCASE 

PRINT INPUT LOAD CASE (TENSION, BENDING, PIN, WEDGE): T, B, P or W> 

READ 1, LCASE 
1 FORMAT (Al) 

PRINT *, ’ INPUT: r/t , b/t , z/t, r/w ?' 

READ *, RT , BT , ZT , RW 

CALL SCF3D ( RT , BT , ZT , RW , LCASE , SCF ) 

PRINT * , ' r/t = ’ ,RT, ’ b/t = J , BT , , z/t = , , ZT , , r/w = ’ ,RW 

PRINT *, Stress-Concentration Factor = SCF 

STOP 

END 

SUBROUTINE SCF3D (RT , BT , ZT , RW , LCASE , SCF) 

C 

C 

C SCF3D - VERSION CREATED APRIL 1991 
C 

C DEVELOPED BY: Kunigal N. Shivakumar and J. C. Newman, Jr. 

C 

C THREE-DIMENSIONAL STRESS-CONCENTRATION IN COUNTERSUNK AND STRAIGHT 
C SHANK RIVET HOLES. 

C 

C PARAMETERS IN THE CALL STATEMENT: 

C (A) INPUT 

C 

C BT = b/t, STRAIGHT SHANK LENGTH TO PLATE THICKNESS RATIO 

C RT = r/t, HOLE RADIUS TO THICKNESS RATIO 

C ZT = z/t, LOCATION WHERE STRESS-CONCENTRATION FACTOR IS REQUIRED 

C AS A RATIO OF PLATE THICKNESS 

C NOTE: IS MEASURED FROM THE MID-PLANE OF THE PLATE 

C RW = r/w, HOLE RADIUS TO PLATE WIDTH 

C LCASE - LOADING CASE AS DEFINED BELOW 

C 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


STRAIGHT SHANK 

LCASE = T, REMOTE TENSION 
LCASE = B, REMOTE BENDING 
LCASE = P, PIN LOADING 
LCASE = W, WEDGE LOADING 

COUNTER SUNK HOLE 

LCASE = T, REMOTE TENSION 
LCASE = B, REMOTE BENDING 

(B) OUTPUT 

SCF - THREE-DIMENSIONAL STRESS-CONCENTRATION FACTOR 

REMOTE TENSION: SCF = MAX. STRESS/S 

S = REMOTE APPLIED STRESS 

REMOTE BENDING: SCF = MAX. STRESS/ (6M/ (t*t) ) 

M = REMOTE APPLIED MOMENT PER 
UNIT WIDTH 

PIN LOAD: SCF = MAX. STRESS/ (P/ (2rt) ) 

P = PIN LOAD 


C WEDGE LOAD: SCF = MAX. STRESS/ (P/ (2rt) ) 

C P = WEDGE LOAD 

C 

C NOTE: SCF FOR SIMULATED RIVET LOADING IS OBTAINED BY 

C ADDING ONE-HALF OF THE SCF FOR REMOTE TENSION 

C S = P/(2wt) AND ONE-HALF OF THE SCF FOR WEDGE 

C LOADING (2w IS TOTAL WIDTH OF PLATE) . 

c 

CHARACTER *1, LCASE, MCASE,NCASE 
DATA MCASE , NCASE/ 1HT , 1HW/ 

C 

C BEGIN ANALYSIS 
C 


C INPUT ERROR WHEN IERR .NE. 0 (PARAMETER OUT OF RANGE) 

I ERR = 0 

IF (BT .LT. 0.0 .OR. BT .GT. 1.0) IERR = 1 
IF (IERR .EQ. 1) PRINT *,*INPUT PARAMETER b/t OUT OF RANGE* 
IF (RT .GT. 2.5) IERR = 2 
IF (RT .LT. .25) IERR = 2 

IF (IERR .EQ. 2) PRINT *,*INPUT PARAMETER r/t OUT OF RANGE* 
IF(ZT .LT. -0.501 .OR. ZT .GT. 0.501) IERR = 3 
IF (IERR .EQ. 3) PRINT *,*INPUT PARAMETER z/t OUT OF RANGE * 
IF (RW .GT. 0.25) IERR = 4 
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IF(IERR .EQ. 4) PRINT *, ’INPUT PARAMETER r/w OUT OF RANGE’ 
IF(IERR .NE. 0) STOP 

IF (BT .EQ. 1.0) THEN 

IF (LCASE .EQ. ’P’ .OR. LCASE . EQ. ’p’) GOTO 10 
CALL SSHANK (RT , ZT , LCASE , SCF ) 

GOTO 20 
10 CONTINUE 

CALL SSHANK (RT , ZT , MCASE , SCFT) 

CALL SSHANK (RT,ZT,NCASE,SCFW) 

SCF = (SCFW + RW * SCFT) * 0.5 
20 RETURN 
ELSE 

IF (BT .EQ. 0.0) THEN 

CALL CSHANK(RT,BT,ZT, LCASE, SCF) 

RETURN 

ENDIF 

IF (BT .GT. 0.0 .AND. BT .LE. 0.25) THEN 
CALL NZT(BT,ZT, 0.0,0. 25, ZT1.ZT2) 

CALL CSHANK (RT , 0 . 0 , ZT1 , LCASE , SCF1 ) 

CALL CSHANK (RT , 0 . 25 , ZT2 , LCASE , SCF2) 

SCF = SCF1 + (SCF2-SCFD/0.25 * BT 

RETURN 

ENDIF 

IF(BT .GT. 0.25 .AND. BT .LE. 0.5) THEN 
CALL NZT(BT,ZT, 0.25,0. 5, ZT1.ZT2) 

CALL CSHANK (RT , 0 . 25 , ZT1 , LCASE , SCF 1 ) 

CALL CSHANK (RT , 0 . 50 , ZT2 , LCASE , SCF2) 

SCF = SCF1 + (SCF2-SCF1 ) /O . 25 * (BT-0.25) 

RETURN 

ENDIF 

IF(BT .GT. 0.50 .AND. BT .LE. 0.75) THEN 
CALL NZT(BT,ZT, 0.5,0. 75, ZT1.ZT2) 

CALL CSHANK(RT,0.50,ZT1,LCASE,SCF1) 

CALL CSHANK (RT , 0 . 75 , ZT2 , LCASE , SCF2) 

SCF = SCF1 + (SCF2-SCFl)/0.25 * (BT-0.50) 

RETURN 

ENDIF 

IF (BT .GT. 0.75 .AND. BT .LT. 1.0) THEN 
CALL NZT (BT, ZT, 0 .75,1.0, ZT1 ,ZT2) 

CALL CSHANK (RT , 0 . 75 , ZT 1 , LCASE , SCF 1 ) 

CALL SSHANK ( RT , ZT , LCASE , SCF2) 

SCF = SCF1 + (SCF2-SCF1) /O . 25 * (BT-0.75) 

RETURN 

ENDIF 



ENDIF 

END 

SUBROUTINE NZT (BT , ZT , BT1 , BT2 , ZT1 , ZT2) 

C 

C EVALUATE APPROPRIATE Z-LOCATION FOR COUNTER-SUNK HOLE 
C 

IF (ZT .GT. (BT-0.5) .AND. ZT .LE. 0.5) THEN 
ZT1 = BT1-0.5 + (ZT-BT+0.5) * (1 . -BT1) / (1 . -BT) 

ZT2 = BT2-0.5 + (ZT-BT+0.5) * (1 . -BT2) / (1 . -BT) 

ELSE 

ZT1 = BT1-0.5 + (ZT-BT+0.5) * BT1/BT 

ZT2 = BT2-0.5 + (ZT-BT+0.5) * BT2/BT 

ENDIF 

RETURN 

END 

SUBROUTINE SSHANK (RT , ZT , LCASE , SCF) 

c 

C 

C THREE-DIMENSIONAL STRESS-CONCENTRATION EQUATION FOR STRAIGHT 
C SHANK RIVET HOLE SUBJECTED TO: 

C (1) REMOTE TENSION 

C (2) REMOTE BENDING 

C (3) PIN LOADING IN HOLE (r/w < 0.25) 

C (4) WEDGE LOADING IN HOLE 

C 

C RANGE OF PARAMETERS: -0.5 < z/t < 0.5; 0.25 < r/t < 2.5 
C 

c 

DIMENSION ALP(4 ,5,2) , ALPB(4,4) 

CHARACTER *1, LCASE 


DATA ALP/ 


1 

3.1825, 

. 1679, 

-.2063, 

.0518, 

2 

.4096, 

-1.5125, 

1.1650, 

-.2539, 

3 

-1.2831, 

2.8632, 

-2.0000, 

.4239, 

4 

2.2778, 

-6.0148, 

4.5357, 

-.9983, 

5 

-2.0712, 

5.2088, 

-3.8337, 

.8331, 

6 

1.7130, 

. 1390, 

-.1356, 

.0317, 

7 

.3626, 

-1.0206, 

.7242, 

-.1527, 

8 

-1.5767, 

3.0242, 

-2.0075, 

.4169, 

9 

3.1870, 

-6.5555, 

4.4847, 

-.9450, 

C 

-2.3673, 

4.6981, 

-3.1644, 

.6614/ 

DATA ALPB/ 




1 

3.1773, 

-1.7469, 

.9801, 

-.1875, 

2 

-.2924, 

.1503, 

-.0395, 

.0040, 

3 

.8610, 

-2.1651, 

1.5684, 

-.3370, 

4 

-1.2427, 

2.7202, 

-1.8804, 

.3957/ 
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IF (LCASE 

.EQ 

j *p > 

.OR. 

LCASE . 

EQ. 

’f) 

i 

ii 

-j 

IF (LCASE 

.EQ 

. 'W 1 

.OR. 

LCASE . 

EQ. 

’w>) 

L = 2 

IF (LCASE 

.EQ 

. *B * 

.OR. 

LCASE . 

EQ. 

’b’) 

CO 

ii 

-J 

IF (L .LT. 

1 . 

OR. L 

.GT. 

3) PRINT *, 

1 LOAD TYPE NOT DEFINED 

IF(L .LT. 

1 . 

OR. L 

.GT. 

3) STOP 




SCF = 0.0 
Z2T = 2 * ZT 
IF(L.LE.2) THEN 
DO 11 I = 1, 4 
II = 1-1 

IF (Z2T .EQ. 0.0) THEN 

SCF = SCF + ALP ( I , 1 , L) * RT**I1 

ELSE 

DO 10 J = 1, 5 

J1 = ( J-l) *2 

SCF = SCF + ALP (I , J , L) * RT**I1 * Z2T**J1 

10 CONTINUE 
ENDIF 

11 CONTINUE 
ELSE 

DO 20 I = 1, 4 
II = I - 1 
DO 20 J = 1, 4 
J1 = 2*J - 1 

SCF = SCF + ALPB (I , J) * RT**I1 * Z2T**J1 
20 CONTINUE 

ENDIF 
RETURN 
END 

SUBROUTINE CSHANK (RT , BT , ZT , LCASE , SCF) 

C 

C 

C THREE-DIMENSIONAL STRESS CONCENTRATION FACTOR FOR COUNTER-SUNK 
C RIVET HOLE SUBJECTED TO: 

C (1) REMOTE TENSION 

C (2) REMOTE BENDING 

C SOLUTION IS FOR THE COUNTER-SUNK ANGLE OF 100 DEGREES AND A 

C SELECTED VALUE OF b/t (RATIOS ARE 0, .25, .50 AND 0.75). 

C RESULTS FOR ANY OTHER b/t VALUE ARE COMPUTED BY LINEAR 

C INTERPOLATION BETWEEN THE TWO LIMITING b/t VALUES. 

C 

C RANGE OF PARAMETERS: -0 . 5 < z/t < 0 . 5 ; 0.25 < r/t < 2.5 
C 

c 

C 

DIMENSION ALP (3 ,5,4,2) , BET(3,5,4,2) 

CHARACTER *1, LCASE 



c 

3.1675, 

1.2562, 

-.4052, 

3.7503, 

-8.8507, 

2.8948,- 

15 . 6036 , 

c 

23.4071, 

-7.7898, 

22.1981,- 

-30.9691, 

10.3670, 

-11.1465, 

15.1933 

c 

-5.0730, 

3.5507, 

.7198, 

-.2232, 

. 1185, 

1.0574, 

-.2623, 

c 

-2.2035, 

2.0077, 

-.4746, 

-4.2715, 

5.0031, 

-1.4629, 

-2.9410 

c 

3.7985, 

-1.1888, 

3.4454, 

.4835, 

-.1485, 

.3460, 

.1089, 

c 

.0844, 

-2.2150, 

1.1287, 

-.0843, 

-6.5876, 

7.3731, 

-2 . 1234. 

c 

-4.9136, 

6.1237, 

-1.8862, 

3.3341, 

.0777, 

-.0259, 

- . 0229 . 

c 

-.5498, 

.3049, 

-4.7184, 

2.8236, 

-.5229, 

-12.1049, 

12.3213 

c 

-3.5036, 

-8.1604, 

9.1806, 

-2.7318, 




c 

-2.7192, 

.4773, 

-.1620, 

5.2713, 

2.1888, 

-.6093, 

3.2839, 

c - 

-10.8632, 

3.2768, 

-4.5453, 

9.9384, 

-3.0428, 

.7327, 

-2.2565 

c 

.7056, 

-1.4221, 

.4322, 

-.1424, 

1.6817, 

-1.1265, 

.3481, 

c 

1.2863, 

-2.0711, 

.6784, 

2.4568, 

-3.8178, 

1.2723, 

1 . 8492 , 

c 

-2.6911, 

.8911, 

.1935, 

-.0883, 

.0135, 

3 . 8939 , 

-2.7731, 

c 

.8887, 

3.2128, 

-6.4904, 

2.3056, 

5.8885, 

-10.6559, 

3. 7525 : 

c 

3.9311, 

-6.2356, 

2.1384, 

1.7020, 

-.7146, 

.2021, 

6.4706, 

c 

-4.6850, 

1.4482, 

8.3737,- 

-12.5101, 

4.0720, 

14.4058,- 

-19.9993 

c 

6.4740, 

8.3649, 

-10.8222, 

3.4552/ 




DATA BET/ 







c 

3.1675, 

1.2562, 

-.4052, 

3.7503, 

-8.8507, 

2.8948,- 

-15.6036 

c 

23.4071, 

-7.7898, 

22.1981,- 

-30.9691, 

10.3670, 

-11 . 1465, 

15.1933 

c 

-5.0730, 

3.5507, 

.7198, 

-.2232, 

-1.4878, 

-4.1557, 

1.2616, 

c 

.6958, 

8.9708, 

-2.6866, 

2.6002, 

-13.8774, 

4.2240, 

-3.0363 

c 

8.2145, 

-2.5264, 

3.4454, 

.4835, 

-.1485, 

-1.1969, 

-2.6156 

c 

.7803, 

1.0127, 

1.8286, 

-.5102, 

.3438, 

-1.8037, 

.5698, 

c 

-1.3109, 

1.7708, 

-.5768, 

3.3341, 

.0777, 

- .0259, 

-.6655, 

c 

-1.7805, 

.5880, 

-.9018, 

3.0805, 

-1.0493, 

2.1386, 

-4.3757 

c 

1.5303, 

-1.6774, 

2.7382, 

-.9445, 




c 

-2.7192, 

.4773, 

-.1620, 

5.2713, 

2.1888, 

- .6093, 

3.2839, 

c - 

-10.8632, 

3.2768, 

-4.5453, 

9.9384, 

-3.0428, 

.7327, 

-2.2565 

c 

.7056, 

-1.4221, 

.4322, 

-.1424, 

6.6870, 

-2.1064, 

.7330, 

c 

-9.2419, 

4.3538, 

-1.5784, 

13.6204, 

-9.2163, 

3.1486, 

-7.6364 

c 

6.0611, 

-2.0053, 

.1935, 

-.0883, 

.0135, 

2.8201, 

-1.4920, 

c 

.5510, 

-.4453, 

1.8097, 

-.7420, 

.6186, 

-3.4144, 

1.2552, 

c 

-1.1330, 

2.6470, 

-.8987, 

1.7020, 

-.7146, 

.2021, 

.2472, 

c 

-.4422, 

.2356, 

1.8402, 

.0875, 

-.2380, 

-1.9081, 

-.4494, 

c 

.4036, 

.1992, 

.8738, 

-.3866/ 




*** BEGIN ANALYSIS 






IF 

(LCASE . 

EQ. . 

OR. LCASE 

.EQ. J t> 

) L = 1 



IF 

(LCASE . 

EQ. . 

OR. LCASE 

. EQ . 5 b ’ 

) L = 2 




IF (L .LT. 1 -OR. L .GT. 2) PRINT *,’LOAD TYPE NOT DEFINED’ 
IF (L .LT. 1 .OR. L . GT. 2) STOP 



IF (LCASE .EQ. ’W’) PRINT *, ’WEDGE LOAD SOLUTION NOT AVAILABLE’ 

IF (LCASE .EQ. ’W’) STOP 

IF (BT .EQ. O.O) K = 1 

IF(BT .EQ. 0.25) K = 2 

IF (BT .EQ. 0.50) K = 3 

IF (BT .EQ. 0.75) K = 4 

CCOR = (1. - 2*BT) /2 . 

SCF = 0.0 

IF (BT .NE. 0.0) THEN 

IF (-0.5 .LE. ZT .AND. ZT .LE. -CCOR) THEN 
T1 = CCOR/BT 
Z = T1 + ZT/BT 
DO 11 I = 1, 3 
II = 1-1 

IF ( Z .EQ. 0.0) THEN 

SCF = SCF + ALP (I , 1 , K , L) * RT**I1 

ELSE 

DO 10 J = 1, 5 
J1 = J - 1 

SCF = SCF + ALP (I , J , K , L) * RT**I1 * Z**J1 

10 CONTINUE 
ENDIF 

11 CONTINUE 
ELSE 

T2 = CCOR/ ( 1 . -BT) 

Z = T2 + ZT/ ( 1 . -BT) 

DO 20 I = 1, 3 
II = 1-1 

IF (Z .EQ. 0.0) THEN 

SCF = SCF + BET (I , 1 ,K,L) * RT**I1 

ELSE 

DO 21 J = 1, 5 
J1 = J - 1 

SCF = SCF + BET (I , J ,K,L) * RT**I1 * Z**J1 
21 CONTINUE 
ENDIF 

20 CONTINUE 
ENDIF 
ELSE 

T2 = CCOR/ (1 . -BT) 

Z = T2 + ZT/ ( 1 . -BT) 

DO 30 I - 1, 3 
II = 1-1 

IF (Z .EQ. 0.0) THEN 

SCF = SCF + BET(I , 1 , K ,L) * RT**I1 

ELSE 
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DO 31 J = 1, 5 
J1 = J - 1 

SCF = SCF + BET (I , J , K , L) * RT**I1 * Z**J1 
31 CONTINUE 
END IF 

30 CONTINUE 
ENDIF 
RETURN 

END 
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Table 1. Details of Finite-Element Idealization of Plates With Countersunk Holes 


b/t 

Nodes 

Elements 

Number 
of layers 

Layer of thickness as percentage 
of plate thickness" 

0 

5756 

1170 

9 

2. 5. 15. 15, 26, 15, 15, 5, 2 

.25 

6938 

1430 

11 

2. 5, 11. 5, 2, 2, 5, 15, 40, 9, 4 

.50 

7529 

1560 

12 

3. 8. 17. 15, 5, 2, 2, 5, 15, 17, 8, 3 

.75 

6938 

1430 

11 

4. 9, 40, 15, 5, 2, 2, 5, 11, 5, 2 


w Layers are numbered from z/t — —0.5 plane. 


Table 2. Comparison of Stress-Concentration Fac tors for Countersunk Holes in Plates 

[b/t = 0.6; 0 r = 90° ] 







Stress-concentration factor from 

Load 

r ft 

w/r 

tfr 

z/t 

Cheng (ref. 17) 

Present 

Tension 
(model 7) 

0.24 

7.2 

7.2 

0.1 

4.19 

"4.310 

Bending 

0.19 

8.2 

8.2 

-0.50 

-2.61 

k— 2.510 

(model 8) 

.19 

8.2 

8.2 

.10 

1.08 

2.216 


.19 

8.2 

8.2 

.50 

2.00 

1.689 


a K tnm = 4.37 at z/t = 0.08. 
h K b = -2.55 at z/t = -0.18. 

"max ' 
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Table 3. Coefficients of Kt Equation for Straight-Shank Hole in Plate 
Subjected to Remote Tension 



Qjj for j of 

i 

0 

2 

4 

6 

8 

0 

3.1825 

0.4096 

-1.2831 

2.2778 

-2.0712 

1 

.1679 

-1.1525 

2.8632 

— 6.0148 

5.2088 

2 

-.2063 

1.1650 

-2.0000 

4.5357 

-3.8337 

3 

.0518 , 

-.2539 

.4239 

-.9983 

.8331 


Table 4. Coefficients of K w Equation for Straight-Shank Hole in Plate 
Subjected to Wedge Loading 



<*ij for j of— 

i 

0 

2 

4 

6 

8 

0 

1.7130 

0.3626 

-1.5767 

3.1870 

-2.3673 

1 

.1390 

-1.0206 

3.0242 

-6.5555 

4.6981 

2 

-.1356 

.7242 

-2.0075 

4.4847 

-3.1644 

3 

.0317 

-.1527 

.4169 

-.9450 

.6614 


Table 5. Coefficients of K ^ Equation for Straight-Shank Hole in Plate 
Subjected to Remote Bending 


i 


<*ij for 

j of 


1 

3 

5 

7 

0 

3.1773 

0.2924 

0.8610 

-1.2427 

1 

-1.7469 

.1503 

-2.1651 

2.7202 

2 

.9801 

-.0395 

1.5684 

-1.8804 

3 

-.1875 

.0040 

-.3370 

.3957 






Table 6. Coefficients of Kf Equation for Countersunk Hole in Plate Subjected to Remote Tension Stress S 


a tJ 

for -0.5 < 

z/t < -(0.5 

- b/t) and j of 

0 

1 

2 

3 

4 

3.1675 

1.2562 

-.4052 

3.7503 

-8.8507 

2.8948 

-15.6036 

23.4071 

-7.7898 

22.1981 

-30.9691 

10.3670 

-11.1465 

15.1933 

-5.0730 

3.5507 

.7198 

-.2232 

0.1185 

1.0574 

-.2623 

-2.2035 
2.0077 
- .4746 

-4.2715 

5.0031 

-1.4629 

-2.9410 

3.7985 

-1.1888 

3.4454 

.4835 

-.1485 

0.3460 

.1089 

.0844 

-2.2150 ' 
1.1287 
-.0843 

-6.5876 

7.3731 

-2.1234 

-4.9136 

6.1237 

-1.8862 

3.3341 

.0777 

-.0259 

-0.0029 

-.5498 

.3049 

-4.7184 
2.8236 
- .5229 

-12.1049 

12.3213 

-3.5036 

-8.1604 

9.1806 

-2.7318 


(0.5 — b/t ) < z/t < 0.5 and j of 


3.1675 
1.2562 
-.4052 
3.5507 
.7108 
-.2232 
3.4454 
.4835 
-.1485 
3.334 f 
.0777 
-.0250 


3.7503 
8.8507 
2.8948 
U ,1878 
1.1557 
1.2016 
hTTogo" 

-2.0156 

.7H03_ 

0.6655” 

-1.7805 

.5880 


-15.6036 
23.4071 
-7.7898 
0.6958 
8.9708 
- 2.6866 
1.0127 
1 .8286 
-.5102 
-0.9018 
3.0805 
-1.0493 


3 

22 1981 
-30.9691 
10.3670 
2.6002 
-13.8774 
4.2240 
0.3438 
-1.8037 
.5698 
2.1386 
-4.3757 
1.5303 


-11.1465 

15.1933 

-5.0730 

-3.0363 

8.2145 

-2.5264 

-1.3109 

1.7708 

-.5768 

-1.6774 

2.7382 

-.9445 


Table 7. Coefficients of Kf, Equation for Countersunk Hole in Plate Subjected to Unit Bending Stress 


2 

050 0~ 


a tJ 

for -0.5 < 

z/t < -(0.5 

- b/t) and j of 

■c f 

7>r —(0.5 — 

b/t) < z/t < 0.5 and j 0 

f 

0 

1 

2 

3 

4 

0 

1 

2 

3 

4 

-2.7192 

5.2713 

3.2839 

-4.5453 

0.7327 

-2.7192 

5.2713 

3.2839 

-4.5453 

0.7327 

.4773 

2.1888 

-10.8632 

9.9384 

-2.2565 

.1773 

2.1888 

-10.8632 

9.9384 

- 2.2565 

-.1620 

-.6093 

3.2768 

, -3.0428 

i 

.7056 

-.1620 

.6093 

3.2768 

-3.0428 

.7056 

-1.4221 

1.6817 

1.2863 

2.4568 

1 .8492 

-1,4221 

6.6870 

-9.2419 

13.6204 

-7.6364 

.4322 

-1.1265 

-2.0711 

-3.8178 

-2.6911 

,4322 

-2.1064 

4.3538 

-9.2163 

6.0611 

-.1424 

.3481 

.6784 

1.2723 

.8911 

-.1424 

.7330 

-1.5784 

3.1486 

-2.0053 

0.1935 

3.8939 

3.2128 

5.8885 

3.9311 

0.1935 

2.8201 

-0.4453 

0.6186 ’ 

-1.1330 

-.0883 

-2.7731 

-6.4904 

-10.6559 

-6.2356 

-.0883 

1,4920 

1.8097 

-3.4144 

2.6470 

.0135 

.8887 

2.3056 

3.7525 

2.1384 

.0135 

.5510 

-.7420 

1.2552 

-.8987 

1.7020 

6.4706 

8.3737 

14.4058 

8.3649 

1.7020 

0.2472 

1.8402 

-1.9081 

0.1992 

-.7146 

-4.6850 

-12.5101 

-19.9993 

-10.8222 

-.7146 

-.4422 

.0875 

-.4494 

.8738 

.2021 

1.4482 

4.0720 

6.4740 

3.4552 

.2021 

.2356 

-.2380 

.4036 

-.3866 





s 

Tension, S 


Bending, M 



(a) Loading oil straight-shank hole. 




(b) Loading on countersunk hole. 

Figure 2. Types of loading on straight-shank and countersunk holes in plates. 



(a) Straight-shank hole; r/t, = 1. 



(b) Countersunk hole; r/t = 0.25; b/t = 0.50; 0 C = ]00°. 

Figure 3. Typical finite-element models for plates with straight-shank and countersunk holes. 
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Figure 10. Comparison of bending SCF equation results with F-E results for straight-shank holes. 




Figure 11. Three-dimensional SCF for pin-loaded plates with straight-shank holes. 













P/2wt P/2wt 

(a) Pin loading. (b) Pin loading. (c) Wedge loading. (d) Tension. 

Figure 14. Superposition of wedge loading and remote tension solutions to calculate pin-load SCF. 
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